Maternal iron status in early pregnancy and DNA methylation in offspring: an epigenome-wide meta-analysis.

BACKGROUND: Unbalanced iron homeostasis in pregnancy is associated with an increased risk of adverse birth and childhood health outcomes. DNA methylation has been suggested as a potential underlying mechanism linking environmental exposures such as micronutrient status during pregnancy with offspring health. We performed a meta-analysis on the association of maternal early-pregnancy serum ferritin concentrations, as a marker of body iron stores, and cord blood DNA methylation. We included 1286 mother-newborn pairs from two population-based prospective cohorts. Serum ferritin concentrations were measured in early pregnancy. DNA methylation was measured with the Infinium HumanMethylation450 BeadChip (Illumina). We examined epigenome-wide associations of maternal early-pregnancy serum ferritin and cord blood DNA methylation using robust linear regression analyses, with adjustment for confounders and performed fixed-effects meta-analyses. We additionally examined whether associations of any CpGs identified in cord blood persisted in the peripheral blood of older children and explored associations with other markers of maternal iron status. We also examined whether similar findings were present in the association of cord blood serum ferritin concentrations with cord blood DNA methylation. RESULTS: Maternal early-pregnancy serum ferritin concentrations were inversely associated with DNA methylation at two CpGs (cg02806645 and cg06322988) in PRR23A and one CpG (cg04468817) in PRSS22. Associations at two of these CpG sites persisted at each of the follow-up time points in childhood. Cord blood serum ferritin concentrations were not associated with cord blood DNA methylation levels at the three identified CpGs. CONCLUSION: Maternal early-pregnancy serum ferritin concentrations were associated with lower cord blood DNA methylation levels at three CpGs and these associations partly persisted in older children. Further studies are needed to uncover the role of these CpGs in the underlying mechanisms of the associations of maternal iron status and offspring health outcomes.

Maternal depression during pregnancy and cord blood DNA methylation: findings from the Avon Longitudinal Study of Parents and Children.

Up to 13% of women may experience symptoms of depression during pregnancy or in the postpartum period. Depression during pregnancy has been associated with an increased risk of adverse neurodevelopmental outcomes in the child and epigenetic mechanisms could be one of the biological pathways to explain this association. In 844 mother-child pairs from the Avon Longitudinal Study of Parents and Children, we carried out an epigenome-wide association study (EWAS) to investigate associations between prospectively collected data on maternal depression ascertained by the Edinburgh Postnatal Depression Scale in pregnancy and DNA methylation in the cord blood of newborn offspring. In individual site analysis, we identified two CpG sites associated with maternal depression in the middle part of pregnancy. In our regional analysis, we identified 39 differentially methylated regions (DMRs). Seven DMRs were associated with depression at any time point during pregnancy, 7 associated with depression in mid-pregnancy, 23 were associated with depression in late pregnancy, and 2 DMRs were associated with depression throughout pregnancy. Several of these map to genes associated with psychiatric disease and brain development. We attempted replication in The Generation R Study and could not replicate our results. Although our findings in ALSPAC suggest that maternal depression could be associated with cord blood DNA methylation the results should be viewed as preliminary and hypothesis generating until further replicated in a larger sample.

Proton therapy of prostate cancer by anterior-oblique beams: implications of setup and anatomy variations.

Proton therapy of prostate by anterior beams could offer an attractive option for treating patients with hip prosthesis and limiting the high-dose exposure to the rectum. We investigated the impact of setup and anatomy variations on the anterior-oblique (AO) proton plan dose, and strategies to manage these effects via range verification and adaptive delivery. Ten patients treated by bilateral (BL) passive-scattering proton therapy (79.2 Gy in 44 fractions) who underwent weekly verification CT scans were selected. Plans with AO beams were additionally created. To isolate the effect of daily variations, initial AO plans did not include range uncertainty margins. The use of fixed planning margins and adaptive range adjustments to manage these effects was investigated. For each case, the planned dose was recalculated on weekly CTs, and accumulated on the simulation CT using deformable registration to approximate the delivered dose. Planned and accumulated doses were compared for each scenario to quantify dose deviations induced by variations. The possibility of estimating the necessary range adjustments before each treatment was explored by simulating the procedure of a diode-based in vivo range verification technique, which would potentially be used clinically. The average planned rectum, penile bulb and femoral heads mean doses were smaller for initial AO compared to BL plans (by 8.3, 16.1 and 25.9 Gy, respectively). After considering interfractional variations in AO plans, the target coverage was substantially reduced. The maximum reduction of V 79.2/D 95/D mean/EUD for AO (without distal margins) (25.3%/10.7/1.6/4.9 Gy, respectively) was considerably larger than BL plans. The loss of coverage was mainly related to changes in water equivalent path length of the prostate after fiducial-based setup, caused by discrepancies in patient anterior surface and bony-anatomy alignment. Target coverage was recovered partially when using fixed planning margins, and fully when applying adaptive range adjustments. The accumulated organs-at-risk dose for AO beams after range adjustment demonstrated full sparing of femoral heads and superior sparing of penile bulb and rectum compared to the conventional BL cases. Our study indicates that using AO beams makes prostate treatment more susceptible to target underdose induced by interfractional variations. Adaptive range verification/adjustment may facilitate the use of anterior beam approaches, and ensure adequate target coverage in every fraction of the treatment.

Stopping Criteria for Log-Domain Diffeomorphic Demons Registration: An Experimental Survey for Radiotherapy Application.

A crucial issue in deformable image registration is achieving a robust registration algorithm at a reasonable computational cost. Given the iterative nature of the optimization procedure an algorithm must automatically detect convergence, and stop the iterative process when most appropriate. This paper ranks the performances of three stopping criteria and six stopping value computation strategies for a Log-Domain Demons Deformable registration method simulating both a coarse and a fine registration. The analyzed stopping criteria are: (a) velocity field update magnitude, (b) mean squared error, and (c) harmonic energy. Each stoping condition is formulated so that the user defines a threshold ∊, which quantifies the residual error that is acceptable for the particular problem and calculation strategy. In this work, we did not aim at assigning a value to e, but to give insights in how to evaluate and to set the threshold on a given exit strategy in a very popular registration scheme. Experiments on phantom and patient data demonstrate that comparing the optimization metric minimum over the most recent three iterations with the minimum over the fourth to sixth most recent iterations can be an appropriate algorithm stopping strategy. The harmonic energy was found to provide best trade-off between robustness and speed of convergence for the analyzed registration method at coarse registration, but was outperformed by mean squared error when all the original pixel information is used. This suggests the need of developing mathematically sound new convergence criteria in which both image and vector field information could be used to detect the actual convergence, which could be especially useful when considering multi-resolution registrations. Further work should be also dedicated to study same strategies performances in other deformable registration methods and body districts.

Radiological Protection in Cone Beam Computed Tomography (CBCT). ICRP Publication 129.

The objective of this publication is to provide guidance on radiological protection in the new technology of cone beam computed tomography (CBCT). Publications 87 and 102 dealt with patient dose management in computed tomography (CT) and multi-detector CT. The new applications of CBCT and the associated radiological protection issues are substantially different from those of conventional CT. The perception that CBCT involves lower doses was only true in initial applications. CBCT is now used widely by specialists who have little or no training in radiological protection. This publication provides recommendations on radiation dose management directed at different stakeholders, and covers principles of radiological protection, training, and quality assurance aspects. Advice on appropriate use of CBCT needs to be made widely available. Advice on optimisation of protection when using CBCT equipment needs to be strengthened, particularly with respect to the use of newer features of the equipment. Manufacturers should standardise radiation dose displays on CBCT equipment to assist users in optimisation of protection and comparisons of performance. Additional challenges to radiological protection are introduced when CBCT-capable equipment is used for both fluoroscopy and tomography during the same procedure. Standardised methods need to be established for tracking and reporting of patient radiation doses from these procedures. The recommendations provided in this publication may evolve in the future as CBCT equipment and applications evolve. As with previous ICRP publications, the Commission hopes that imaging professionals, medical physicists, and manufacturers will use the guidelines and recommendations provided in this publication for implementation of the Commission's principle of optimisation of protection of patients and medical workers, with the objective of keeping exposures as low as reasonably achievable, taking into account economic and societal factors, and consistent with achieving the necessary medical outcomes.

Validation of a deformable image registration technique for cone beam CT-based dose verification.

PURPOSE: As radiation therapy evolves toward more adaptive techniques, image guidance plays an increasingly important role, not only in patient setup but also in monitoring the delivered dose and adapting the treatment to patient changes. This study aimed to validate a method for evaluation of delivered intensity modulated radiotherapy (IMRT) dose based on multimodal deformable image registration (dir) for prostate treatments. METHODS: A pelvic phantom was scanned with CT and cone-beam computed tomography (CBCT). Both images were digitally deformed using two realistic patient-based deformation fields. The original CT was then registered to the deformed CBCT resulting in a secondary deformed CT. The registration quality was assessed as the ability of the dir method to recover the artificially induced deformations. The primary and secondary deformed CT images as well as vector fields were compared to evaluate the efficacy of the registration method and it's suitability to be used for dose calculation. plastimatch, a free and open source software was used for deformable image registration. A B-spline algorithm with optimized parameters was used to achieve the best registration quality. Geometric image evaluation was performed through voxel-based Hounsfield unit (HU) and vector field comparison. For dosimetric evaluation, IMRT treatment plans were created and optimized on the original CT image and recomputed on the two warped images to be compared. The dose volume histograms were compared for the warped structures that were identical in both warped images. This procedure was repeated for the phantom with full, half full, and empty bladder. RESULTS: The results indicated mean HU differences of up to 120 between registered and ground-truth deformed CT images. However, when the CBCT intensities were calibrated using a region of interest (ROI)-based calibration curve, these differences were reduced by up to 60%. Similarly, the mean differences in average vector field lengths decreased from 10.1 to 2.5 mm when CBCT was calibrated prior to registration. The results showed no dependence on the level of bladder filling. In comparison with the dose calculated on the primary deformed CT, differences in mean dose averaged over all organs were 0.2% and 3.9% for dose calculated on the secondary deformed CT with and without CBCT calibration, respectively, and 0.5% for dose calculated directly on the calibrated CBCT, for the full-bladder scenario. Gamma analysis for the distance to agreement of 2 mm and 2% of prescribed dose indicated a pass rate of 100% for both cases involving calibrated CBCT and on average 86% without CBCT calibration. CONCLUSIONS: Using deformable registration on the planning CT images to evaluate the IMRT dose based on daily CBCTs was found feasible. The proposed method will provide an accurate dose distribution using planning CT and pretreatment CBCT data, avoiding the additional uncertainties introduced by CBCT inhomogeneity and artifacts. This is a necessary initial step toward future image-guided adaptive radiotherapy of the prostate.

Validation of automatic contour propagation for 4D treatment planning using multiple metrics.

The aim of this work is to provide insights into multiple metrics clinical validation of deformable image registration and contour propagation methods in 4D lung radiotherapy planning. The following indices were analyzed and compared: Volume Difference (VD), Dice Similarity Coefficient (DSC), Positive Predictive Value (PPV) and Surface Distances (SD). The analysis was performed on three patient datasets, using as reference a ground-truth volume generated by means of Simultaneous Truth And Performance Level Estimation (STAPLE) algorithm from the outlines of five experts. Significant discrepancies in the quality assessment provided by the different metrics in all the examined cases were found. Metrics sensitivity was more evident in presence of image artifacts and particularly for tubular anatomical structures, such as esophagus or spinal cord. Volume Differences did not account for position and DSC exhibited criticalities due to its intrinsic symmetry (i.e. over- and under-estimation of the reference contours cannot be discriminated) and dependency on the total volume of the structure. PPV analysis showed more robust performance, as each voxel concurs to the classification of the propagation, but was not able to detect inclusion of propagated and ground-truth volumes. Mesh distances could interpret the actual shape of the structures, but might report higher mismatches in case of large local differences in the contour surfaces. According to our study, the combination of VD and SD for the validation of contour propagation algorithms in 4D could provide the necessary failure detection accuracy.

Computer models to study uterine activation at labour.

Improving our understanding of the initiation of labour is a major aim of modern obstetric research, in order to better diagnose and treat pregnant women in which the process occurs abnormally. In particular, increased knowledge will help us identify the mechanisms responsible for preterm labour, the single biggest cause of neonatal morbidity and mortality. Attempts to improve our understanding of the initiation of labour have been restricted by the inaccessibility of gestational tissues to study during pregnancy and at labour, and by the lack of fully informative animal models. However, computer modelling provides an exciting new approach to overcome these restrictions and offers new insights into uterine activation during term and preterm labour. Such models could be used to test hypotheses about drugs to treat or prevent preterm labour. With further development, an effective computer model could be used by healthcare practitioners to develop personalized medicine for patients on a pregnancy-by-pregnancy basis. Very promising work is already underway to build computer models of the physiology of uterine activation and contraction. These models aim to predict changes and patterns in uterine electrical excitation during term labour. There have been far fewer attempts to build computer models of the molecular pathways driving uterine activation and there is certainly scope for further work in this area. The integration of computer models of the physiological and molecular mechanisms that initiate labour will be particularly useful.

Interplay effects in proton scanning for lung: a 4D Monte Carlo study assessing the impact of tumor and beam delivery parameters.

Relative motion between a tumor and a scanning proton beam results in a degradation of the dose distribution (interplay effect). This study investigates the relationship between beam scanning parameters and the interplay effect, with the goal of finding parameters that minimize interplay. 4D Monte Carlo simulations of pencil beam scanning proton therapy treatments were performed using the 4DCT geometry of five lung cancer patients of varying tumor size (50.4-167.1 cc) and motion amplitude (2.9-30.1 mm). Treatments were planned assuming delivery in 35 × 2.5 Gy(RBE) fractions. The spot size, time to change the beam energy (τes), time required for magnet settling (τss), initial breathing phase, spot spacing, scanning direction, scanning speed, beam current and patient breathing period were varied for each of the five patients. Simulations were performed for a single fraction and an approximation of conventional fractionation. For the patients considered, the interplay effect could not be predicted using the superior-inferior motion amplitude alone. Larger spot sizes (σ ~ 9-16 mm) were less susceptible to interplay, giving an equivalent uniform dose (EUD) of 99.0 ± 4.4% (1 standard deviation) in a single fraction compared to 86.1 ± 13.1% for smaller spots (σ ~ 2-4 mm). The smaller spot sizes gave EUD values as low as 65.3% of the prescription dose in a single fraction. Reducing the spot spacing improved the target dose homogeneity. The initial breathing phase can have a significant effect on the interplay, particularly for shorter delivery times. No clear benefit was evident when scanning either parallel or perpendicular to the predominant axis of motion. Longer breathing periods decreased the EUD. In general, longer delivery times led to lower interplay effects. Conventional fractionation showed significant improvement in terms of interplay, giving a EUD of at least 84.7% and 100.0% of the prescription dose for the small and larger spot sizes respectively. The interplay effect is highly patient specific, depending on the motion amplitude, tumor location and the delivery parameters. Large degradations of the dose distribution in a single fraction were observed, but improved significantly using conventional fractionation.

Auto-masked 2D/3D image registration and its validation with clinical cone-beam computed tomography.

Image-guided alignment procedures in radiotherapy aim at minimizing discrepancies between the planned and the real patient setup. For that purpose, we developed a 2D/3D approach which rigidly registers a computed tomography (CT) with two x-rays by maximizing the agreement in pixel intensity between the x-rays and the corresponding reconstructed radiographs from the CT. Moreover, the algorithm selects regions of interest (masks) in the x-rays based on 3D segmentations from the pre-planning stage. For validation, orthogonal x-ray pairs from different viewing directions of 80 pelvic cone-beam CT (CBCT) raw data sets were used. The 2D/3D results were compared to corresponding standard 3D/3D CBCT-to-CT alignments. Outcome over 8400 2D/3D experiments showed that parametric errors in root mean square were <0.18° (rotations) and <0.73 mm (translations), respectively, using rank correlation as intensity metric. This corresponds to a mean target registration error, related to the voxels of the lesser pelvis, of <2 mm in 94.1% of the cases. From the results we conclude that 2D/3D registration based on sequentially acquired orthogonal x-rays of the pelvis is a viable alternative to CBCT-based approaches if rigid alignment on bony anatomy is sufficient, no volumetric intra-interventional data set is required and the expected error range fits the individual treatment prescription.

WE-E-213CD-02: Gaussian Weighted Multi-Atlas Based Segmentation for Head and Neck Radiotherapy Planning.

PURPOSE: To develop a multi-atlas segmentation strategy for IMRT head and neck therapy planning. METHODS: The method was tested on thirty-one head and neck simulation CTs, without demographic or pathology pre-clustering. We compare Fixed Number (FN) and Thresholding (TH) selection (based on normalized mutual information ranking) of the atlases to be included for current patient segmentation. Next step is a pairwise demons Deformable Registration (DR) onto current patient CT. DR was extended to automatically compensate for patient different field of view. Propagated labels are combined according to a Gaussian Weighted (GW) fusion rule, adapted to poor soft tissues contrast. Agreement with manual segmentation was quantified in terms of Dice Similarity Coefficient (DSC). Selection methods, number of atlases used, as well as GW, average and majority voting fusion were discriminated by means of Friedman Test (a=5%). Experimental tuning of the algorithm parameters was performed on five patients, deriving an optimal configuration for each structure. RESULTS: DSC reduction was not significant when ten or more atlases are selected, whereas DSC for single most similar atlas selection is 10% lower in median. DSC of FN selection rule were significantly higher for most structures. Tubular structures may benefit from computing average contour rather than looking at the singular voxel contribution, whereas the best performing strategy for all other structures was GW. When half database is selected, final median DSC were 0.86, 0.80, 0.51, 0.81, 0.69 and 0.79 for mandible, spine, optical nerves, eyes, parotids and brainstem respectively. CONCLUSION: We developed an efficient algorithm for multiatlas based segmentation of planning CT volumes, based on DR and GW. FN selection of database atlases is foreseen to increase computational efficiency. The absence of clinical pre-clustering and specific imaging protocol on database subjects makes the results closer to real clinical application. "Progetto Roberto Rocca" funded by the Fondazione Fratelli Agostino and Enrico Rocca, Italy.

Ranking of stopping criteria for log domain diffeomorphic demons application in clinical radiation therapy.

Deformable Image Registration is a complex optimization algorithm with the goal of modeling a non-rigid transformation between two images. A crucial issue in this field is guaranteeing the user a robust but computationally reasonable algorithm. We rank the performances of four stopping criteria and six stopping value computation strategies for a log domain deformable registration. The stopping criteria we test are: (a) velocity field update magnitude, (b) vector field Jacobian, (c) mean squared error, and (d) harmonic energy. Experiments demonstrate that comparing the metric value over the last three iterations with the metric minimum of between four and six previous iterations is a robust and appropriate strategy. The harmonic energy and vector field update magnitude metrics give the best results in terms of robustness and speed of convergence.

On developing B-spline registration algorithms for multi-core processors.

Spline-based deformable registration methods are quite popular within the medical-imaging community due to their flexibility and robustness. However, they require a large amount of computing time to obtain adequate results. This paper makes two contributions towards accelerating B-spline-based registration. First, we propose a grid-alignment scheme and associated data structures that greatly reduce the complexity of the registration algorithm. Based on this grid-alignment scheme, we then develop highly data parallel designs for B-spline registration within the stream-processing model, suitable for implementation on multi-core processors such as graphics processing units (GPUs). Particular attention is focused on an optimal method for performing analytic gradient computations in a data parallel fashion. CPU and GPU versions are validated for execution time and registration quality. Performance results on large images show that our GPU algorithm achieves a speedup of 15 times over the single-threaded CPU implementation whereas our multi-core CPU algorithm achieves a speedup of 8 times over the single-threaded implementation. The CPU and GPU versions achieve near-identical registration quality in terms of RMS differences between the generated vector fields.

Four-dimensional targeting error analysis in image-guided radiotherapy.

Image-guided therapy (IGT) involves acquisition and processing of biomedical images to actively guide medical interventions. The proliferation of IGT technologies has been particularly significant in image-guided radiotherapy (IGRT), as a way to increase the tumor targeting accuracy. When IGRT is applied to moving tumors, image guidance becomes challenging, as motion leads to increased uncertainty. Different strategies may be applied to mitigate the effects of motion: each technique is related to a different technological effort and complexity in treatment planning and delivery. The objective comparison of different motion mitigation strategies can be achieved by quantifying the residual uncertainties in tumor targeting, to be detected by means of IGRT technologies. Such quantification requires an extension of targeting error theory to a 4D space, where the 3D tumor trajectory as a function of time measured (4D Targeting Error, 4DTE). Accurate 4DTE analysis can be represented by a motion probability density function, describing the statistical fluctuations of tumor trajectory. We illustrate the application of 4DTE analysis through examples, including weekly variations in tumor trajectory as detected by 4DCT, respiratory gating via external surrogates and real-time tumor tracking.

A critical role for TRAIL in resolution of granulomatous experimental autoimmune thyroiditis.

Granulomatous experimental autoimmune thyroiditis (G-EAT) is induced by mouse thyroglobulin (MTG)-sensitized splenocytes activated in vitro with MTG and IL-12. Thyroid lesions reach maximal severity 20 days after cell transfer, and usually resolve or progress to fibrosis by day 60 depending on the extent of thyroid damage at day 20. Our previous studies indicated that neutralization of TNF-alpha or FasL had no effect on G-EAT induction, but neutralization of TNF-alpha promoted, while neutralization of FasL inhibited, G-EAT resolution. TNF-related apoptosis-inducing ligand (TRAIL) is a member of the TNF superfamily. This study was undertaken to define the role of endogenous TRAIL in G-EAT development and/or resolution. Neutralization of endogenous TRAIL had little effect on G-EAT induction, but significantly inhibited G-EAT resolution and increased thyroid fibrosis. This correlated with higher expression of pro-inflammatory cytokines and preferential expression of the pro-apoptotic molecule TRAIL, and anti-apoptotic molecules FLIP and Bcl-xL on inflammatory cells in thyroids of anti-TRAIL-treated recipients. The results suggest that endogenous TRAIL is not required for G-EAT development in recipients, but is critical for G-EAT resolution. Endogenous TRAIL might promote resolution, at least in part, through modulation of the balance between pro- and anti-inflammatory cytokines, and the expression pattern of pro- and anti-apoptotic molecules of thyroid epithelial cells (TECs) and inflammatory cells.

GPU-based streaming architectures for fast cone-beam CT image reconstruction and demons deformable registration.

This paper shows how to significantly accelerate cone-beam CT reconstruction and 3D deformable image registration using the stream-processing model. We describe data-parallel designs for the Feldkamp, Davis and Kress (FDK) reconstruction algorithm, and the demons deformable registration algorithm, suitable for use on a commodity graphics processing unit. The streaming versions of these algorithms are implemented using the Brook programming environment and executed on an NVidia 8800 GPU. Performance results using CT data of a preserved swine lung indicate that the GPU-based implementations of the FDK and demons algorithms achieve a substantial speedup--up to 80 times for FDK and 70 times for demons when compared to an optimized reference implementation on a 2.8 GHz Intel processor. In addition, the accuracy of the GPU-based implementations was found to be excellent. Compared with CPU-based implementations, the RMS differences were less than 0.1 Hounsfield unit for reconstruction and less than 0.1 mm for deformable registration.

Derivation of the tumor position from external respiratory surrogates with periodical updating of the internal/external correlation.

In this work we develop techniques that can derive the tumor position from external respiratory surrogates (abdominal surface motion) through periodically updated internal/external correlation. A simple linear function is used to express the correlation between the tumor and surrogate motion. The function parameters are established during a patient setup session with the tumor and surrogate positions simultaneously measured at a 30 Hz rate. During treatment, the surrogate position, constantly acquired at 30 Hz, is used to derive the tumor position. Occasionally, a pair of radiographic images is acquired to enable the updating of the linear correlation function. Four update methods, two aggressive and two conservative, are investigated: (A1) shift line through the update point; (A2) re-fit line through the update point; (C1) re-fit line with extra weight to the update point; (C2) minimize the distances to the update point and previous line fit point. In the present study of eight lung cancer patients, tumor and external surrogate motion demonstrate a high degree of linear correlation which changes dynamically over time. It was found that occasionally updating the correlation function leads to more accurate predictions than using external surrogates alone. In the case of high imaging rates during treatment (greater than 2 Hz) the aggressive update methods (A1 and A2) are more accurate than the conservative ones (C1 and C2). The opposite is observed in the case of low imaging rates.

Effects of organ motion on IMRT treatments with segments of few monitor units.

Interplay between organ (breathing) motion and leaf motion has been shown in the literature to have a small dosimetric impact for clinical conditions (over a 30 fraction treatment). However, previous studies did not consider the case of treatment beams made up of many few-monitor-unit (MU) segments, where the segment delivery time (1-2 s) is of the order of the breathing period (3-5 s). In this study we assess if breathing compromises the radiotherapy treatment with IMRT segments of low number of MUs. We assess (i) how delivered dose varies, from patient to patient, with the number of MU per segment, (ii) if this delivered dose is identical to the average dose calculated without motion over the path of the motion, and (iii) the impact of the daily variation of the delivered dose as a function of MU per segment. The organ motion was studied along two orthogonal directions, representing the left-right and cranial-caudal directions of organ movement for a patient setup in the supine position. Breathing motion was modeled as sin(x), sin4(x), and sin6(x), based on functions used in the literature to represent organ motion. Measurements were performed with an ionization chamber and films. For a systematic study of motion effects, a MATLAB simulation was written to model organ movement and dose delivery. In the case of a single beam made up of one single segment, the dose delivered to point in a moving target over 30 fractions can vary up to 20% and 10% for segments of 10 MU and 20 MU, respectively. This dose error occurs because the tumor spends most of the time near the edges of the radiation beam. In the case of a single beam made of multiple segments with low MU, we observed 2.4%, 3.3%, and 4.3% differences, respectively, for sin(x), sin4(x), and sin6(x) motion, between delivered dose and motion-averaged dose for points in the penumbra region of the beam and over 30 fractions. In approximately 5-10% of the cases, differences between the motion-averaged dose and the delivered 30-fraction dose could reach 6%, 8% and 10-12%, respectively for sin(x), sin4(x), and sin6(x) motion. To analyze a clinical IMRT beam, two patient plans were randomly selected. For one of the patients, the beams showed a likelihood of up to 25.6% that the delivered dose would deviate from the motion-averaged dose by more than 1%. For the second patient, there was a likelihood of up to 62.8% of delivering a dose that differs by more than 1% from the motion-averaged dose and a likelihood of up to approximately 30% for a 2% dose error. For the entire five-beam IMRT plan, statistical averaging over the beams reduces the overall dose error between the delivered dose and the motion-averaged dose. For both patients there was a likelihood of up to 7.0% and 33.9% that the dose error was greater than 1%, respectively. For one of the patients, there was a 12.6% likelihood of a 2% dose error. Daily intrafraction variation of the delivered dose of more than 10% is non-negligible and can potentially lead to biological effects. We observed [for sin(x), sin4(x), and sin6(x)] that below 10-15 MU leads to large daily variations of the order of 15-35%. Therefore, for small MU segments, non-negligible biological effects can be incurred. We conclude that for most clinical cases the effects may be small because of the use of many beams, it is desirable to avoid low-MU segments when treating moving targets. In addition, dose averaging may not work well for hypo-fractionation, where fewer fractions are used. For hypo-fractionation, PDF modeling of the tumor motion in IMRT optimization may not be adequate.

How much margin reduction is possible through gating or breath hold?

We determined the relationship between intra-fractional breathing motion and safety margins, using daily real-time tumour tracking data of 40 patients (43 tumour locations), treated with radiosurgery at Hokkaido University. We limited our study to the dose-blurring effect of intra-fractional breathing motion, and did not consider differences in positioning accuracy or systematic errors. The additional shift in the prescribed isodose level (e.g. 95 %) was determined by convolving a one-dimensional dose profile, having a dose gradient representing an 8 MV beam through either lung or water, with the probability density function (PDF) of breathing. This additional shift is a measure for the additional margin that should be applied in order to maintain the same probability of tumour control as without intra-fractional breathing. We show that the required safety margin is a nonlinear function of the peak-to-peak breathing motion. Only a small reduction in the shift of isodose curves was observed for breathing motion up to 10 mm. For larger motion, 20 or 30 mm, control of patient breathing during irradiation, using either gating or breath hold, can allow a substantial reduction in safety margins of about 7 or 12 mm depending on the dose gradient prior to blurring. Clinically relevant random setup uncertainties, which also have a blurring effect on the dose distribution, have only a small effect on the margin needed for intra-fractional breathing motion. Because of the one-dimensional nature of our analysis, the resulting margins are mainly applicable in the superior-inferior direction. Most measured breathing PDFs were not consistent with the PDF of a simple parametric curve such as cos4, either because of irregular breathing or base-line shifts. Instead, our analysis shows that breathing motion can be modelled as Gaussian with a standard deviation of about 0.4 times the peak-to-peak breathing motion.